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ABSTRACT 



The study represents a continuation of the development of 
a model output statistics scheme to specify marine fog over 
the open ocean and in coastal waters. Thirty-seven direct 
and derived Fleet Numerical Weather Central model output para- 
meters, monthly climatological fog frequencies, combinations 
of the aforementioned parameters (i.e., interactive parameters) 
and a persistence parameter are used as predictors in a step- 
wise multiple linear regression approach to estimate a predic- 
tand defined as marine fog probability. The predictand is 
categorized in two ways, in one case (FOGCAT I) as smoothed 
probabilities from 0 to 100% as a function of present weather, 
past weather, visibility and low cloud type; and, in another 
case (FOGCAT II) as a limited number of discrete probabilities 
(to include 0 and 100%) derived from present weather, past 
weather and visibility only. This study derives diagnostic 
regression equations only using as a dependent data sample 
over 24,000 North Pacific Ocean (30-60N) surface synoptic ship 
observations at 0000 GMT for June through August 1976 and 1977. 
The predictor parameters contributing most significantly to 
the variance are sensible and evaporative heat fluxes, monthly 
climatological fog frequencies, and meridional wind speed. 
Threat, Heidke skill, and Panof sky -Brier probability scoring 
methods are applied to a selection of the derived equations. 
Predictand variance explained reaches .670, threat/ skill 
scores reach 0.42 and probability scores are as low as 0.28 
using the FOGCAT I categorization scheme for the predictand. 
Equations for June and July appear more stable than those for 
August . 
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I. INTRODUCTION 



Marine fog presents a threat to safe nautical activities 
as well as to low-level aviation over the oceans. Histori- 
cal monetary and human losses associated with United States 
Navy operations, attributable solely or mostly to fog, have 
been documented by Wheeler and Leipper (1974) . These types 
of losses should be significantly reduced as the methods for 
analysis and forecasting of marine fog become more accurate. 
The environmental science community would additionally bene- 
fit from the increased accuracy of fog analyses/ forecasts by 
a commensurate improvement in specifying fog-associated para- 
meters in various analysis models, especially boundary layer 
models . 

Research in the past several years conducted by the De- 
partments of Meteorology and Oceanography, Naval Postgraduate 
School CNPS)_, Monterey, California, has been directed at more 
adequate diagnoses of marine fog and a better representation 
of its climatology, as based on ocean-station-vessel and 
ship-of-opportunity observations CR enar d, Englebretson and 
Daughenbaugh , 1975; Renard, 1976). Additionally, more recent 
research, has attempted to use geostationary weather satellite 
data to diagnose marine fog areas (Ihli and Renard, 1977; 
McNab, 1979). The totality and continuity of such satellite 
observations makes this approach potentially very useful. 
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However, much further testing is needed before operational 
application is feasible. 

In a more statistical vein, an earlier study by Nelson 
(1972) developed regression equations using the "perfect 
prognostic approach" to forecast visibility (and hence fog) 
at ocean station vessels in the North Atlantic Ocean. This 
approach is not unlike the Model Output Statistics (MOS) 
methods used currently by the National Weather Service (NWS) 
to forecast weather parameters of interest over the conti- 
nental U.S., Hawaii and Alaska (Glahn and Lowry, 1972), ex- 
cept the "imperfect prognostic approach" is the present mode. 
Nelson found quite high correlations of visibility with wind 
speed, relative moisture content and various evaporative 
parameters, using observed data. However, when he attempted 
to forecast fog using output parameters from the United 
States Navy Fleet Numerical Weather Central's (FNWC) numeri- 
cal prediction models, his regression equations showed little 
skill. Earlier work, conducted by Schramm (1966), found high 
correlations between observed values of humidity, air tempera- 
ture, wind speed and visibility at sea; however, the regres- 
sion equation developed from these parameters was of little value 
at times of low visibility. 

Additional research into the feasibility of developing 
regression equations to forecast marine fog, based on FNWC ' s 
model output parameters, has been conducted by Van Orman and 
Renard 01977) and Quinn 01978). Both studies developed 
schemes which statistically processed up to 37 direct or 
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derived FNWC model output analysis-time'*" parameters in con- 
junction with surface synoptic ship observations to 
generate multiple linear regression equations yielding as a 

predictand the probability of marine fog occurrence at 
analysis time over the North Pacific Ocean. The Van Orman 
/Renard results, based solely on July 1976 data, showed 
considerable skill over FNWC ' s advective fog model, FTER, 
as applied to the analysis of fog (U.S. Naval Weather Service 
Command, 1975) and climatological fog frequencies developed 
at the NPS OVillms, 1975). Quinn expanded the data base to 
also include June and August 1976 and developed regression 
equations both for the unique probabilistic marine fog pre- 
dictand introduced by Van Orman and Renard as well as a 
modification of that predictand. Quinn further introduced 
interactive predictor parameters into the regression scheme 
by combining the climatological fog frequency parameter with 
the two most significant model output parameters to derive a 
new set of diagnostic regression equations. A clear improve- 
ment over climatology and FNWC's FTER was demonstrated using 
these modified regression equations on dependent data. The 
equations based on the June 1976 data also showed improvement 
over FTER and climatology when used on an independent data 
set from July 1976. These prior studies give clear indication 
of the relative accuracy of the MOS approach in diagnosing 
marine fog, with a successful extension to forecast modes yet 
to be demonstrated. 

‘'’Some parameters are taken from analysis models, other are 
diagnostic parameters from prognostic models. 
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Currently, only FNWC is producing fog analysis/forecasts 
on a large scale over the oceanic areas through their sta- 
tistical probabilistic product, FTER. A climatological 
parameter is not used in their product due to the unavail- 
ability of an accurate climatology at the time of its 
development. A comprehensive uniform climatology of fog 
occurrence over all of the Northern Hemisphere is now avail- 
able from the National Climatic Center (Guttman, 1978) and 
has been incorporated into this study. This inclusion re- 
presents a refinement of the fog parameter from the two 
prior studies. 



14 



II. OBJECTIVES AND APPROACH 



The primary objective of this study was to continue the 
development of a multiple linear regression approach begun 
by Van Orman and Renard (1977) and pursued by Quinn (1978) 
to specify effectively the distribution and likelihood of 
marine fog at analysis time over the middle latitude (30°- 
60°N) region of the North Pacific Ocean during the summer 
season. Part of the primary objective was to evaluate the 
equations' skill in comparison to the climatological fog fre- 
quencies derived by Guttman (1978) and FNWC ' s existing opera- 
tional scheme for forecasting advection fog probability, 

FTER. A second objective was to determine whether the pre- 
dictand categorization scheme FOGCAT I developed by Van Orman 
and Renard or the FOGCAT II scheme developed by Quinn is 
superior in defining the occurrence of fog. Additionally, 
an investigation of the usefulness of persistence as a 
diagnostic prediction parameter was initiated. The approach 
follows that of Quinn's (1978) experiments with the addition 
of June, July and August 1977 data. 

Sets of regression equations were developed for both fog 
categorization schemes for (a) each of the six months June, July 
and August 1976 and 1977; (b) the combined data for both 

Junes, both Julys and both Augusts, (c) the combined data 
for June, July and August 1976, (d) the combined data for 
June, July, and August 1977, and (e) all of the data for the 
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six-month period. Those equations which were based on data 
from both June's, both July's, both August's, and all six 
months were scored using three methods, skill, threat and 
probability . 
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III. DATA 



A . AREA 

The area of study was confined to the North Pacific Ocean 
to maintain continuity and consistency with the previous 
studies of Van Orman and Renard (1977) and Quinn (1978), 
specifically the area north of 30°N and south of 60°N (Figure 
1). A 27x14 grid was superimposed on this region with grid 
points coincident with their counterparts on the same subarea 
of the standard FNWC 63x63 grid (Figure 2). 

B. TIME PERIOD 

Due to the high frequency of marine fog in the summer 
months over the North Pacific Ocean (Willms, 1975; Guttman , 
1978), the additional summer months (June, July and August) 
of 1977 were added to the available set of the same months 
in 1976. At 0000 GMT the entire North Pacific Ocean is in 
daylight with local noon occurring at the international date- 
line. Only data from this time were used since the accuracy 
of visibility reports from ships transitting the area should 
be best during daylight hours. 

C. SYNOPTIC WEATHER REPORTS 

The June-August 1976 and 1977 synoptic weather reports 
used in this study were provided by the Naval Weather Service 
Detachment, physically located with the National Oceanic and 
Atmospheric Administration’s National Climatic Center, 
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Asheville, North Carolina. These data had been edited to 
eliminate duplicate reports. The total number of observa- 
tions received were 4277 for June 1976, 4391 for July 1976, 
4134 for August 1976, 4232 for June 1977, 4198 for July 
1977 and 3156 for August 1977. All observations were ob- 
tained from transitting ships and Ocean Weather Station 
Papa except data obtained from 11 land stations located in 
the study area (Appendix A). As all of these land stations 
are located on islands or immediately adjacent to the coast 
and have relatively low elevations, the effect of topography 
is thought to be minimum. 

D. MODEL OUTPUT PARAMETERS 

FNWC provided the 22 diagnostic model output parameters 
(MOP's) for the time period and area of interest as output 
from several of their numerical analysis and prediction 
models, the Mass Structure Analysis Model, the Primitive 
Equation Prediction Model, the Marine Wind Model, and the 
Spectral Ocean Wave Model (U.S. Naval Weather Service Command, 
1975), Additionally, a set of 15 parameters were derived 
from the primary set and eight interactive parameters were 

developed. A complete listing and description of all these 
parameters is located in Appendix B. 

E. FOG FREQUENCY CLIMATOLOGY 

Previous research used the Naval Postgraduate School K s 
North Pacific Ocean marine fog climatology (Willms, 1975). 

As an updated hemispheric fog frequency climatology (Guttman, 
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1978) is now available, it was used exclusively throughout 
this research. Specific values for fog frequency were inter 
polated to an 89x45 grid in the area of study, which is 
approximately four times as dense as the FNWC grid over the 
region. This finer resolution allowed a more accurate spec! 
fication of the fog frequency climatology. The fog frequency 
climatologies for each month are presented in Figures 3-5. 
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IV . PROCEDURE 



A. PREDICTAND CATEGORIZATION METHODS 

The two methods of predictand (fog) categorization used 
throughout this study are based on four elements of the 
synoptic reports, namely: present weather (ww) , past weather 
(W) , visibility (VIS), and low cloud type (CL). Observation 
elements other than present weather are used due to incon- 
sistencies in the reported observations, which result in 
part from current rules for coding observed data (U . S . De- 
partments of Commerce, Defense and Transportation, 1969) 

(See Appendix C). For example, fog is not coded as present 
weather whenever any form of precipitation is occurring 
simultaneously. Both methods of fog categorization assign 
a probability of fog occurrence to each synoptic report, 
differing mainly by the number and combination of observa- 
tion elements used to assign the probability value. 

The Fog Categorization I (FOGCAT I) scheme was developed 
by Van Orman and Renard (1977) and used all of the four 
elements listed above. This scheme assigns one of five 
major fog categories and one of three subcategories to each 
synoptic report as a function of ww, W, and CL. Then a fog 
probability, 0 to 100% is assigned to each observation de- 
pending upon the major category and subcategory assigned. 
Appendix D outlines FOGCAT I. 
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The Fog Categorization II (FOGCAT II) scheme developed 
by Quinn (1978) uses only present weather, past weather and 
visibility to assign one of seven discrete fog probabilities 
(0%, 10%, 15%, 35%, 60%, 85% and 100%) to each synoptic 
observation. This scheme is more direct in its method of 
categorization and assigns fog probabilities based on a 
filtering technique. First present weather is evaluated. 

If present weather definitely indicates either the presence 
or absence of fog, then a probability of 100% or 0% is 
assigned, respectively. In either of these cases, both past 
weather and visibility are ignored. For all other present 
weather codes, past weather and visibility are considered 
in the assignment of a fog probability. Past weather is 
evaluated first. If the possibility of fog is eliminated by 
the past weather code, a score of 0% is assigned and 
visibility is not considered. If the past weather code in- 
dicates the possibility of fog, then the visibility is taken 
into account in assigning one of the choices for an inter- 
mediate fog probability (i . e . 10, 15, 35, 60, 80). The 
assignment of these intermediate values is somewhat arbitrary. 
Details of FOGCAT II are presented in Appendix E. 

Neither the FOGCAT I nor FOGCAT II schemes represent an 
ideal categorization scheme for fog probability. FOGCAT I, 

on the one hand, could assign an observation a fog probability 
as low as 62.1% even though the present weather code definitely 
indicates the presence of fog, or it could assign a proba- 
bility as high as 24.1% when fog appears very improbable from 
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the synoptic report. FOGCAT II was developed to achieve a 
better fit of probability to occurrence of event. Except 
for the "100" and "0" values, the FOGCAT II probabilities, 
like those of FOGCAT I, are somewhat arbitrary. In any case 
the schemes were accepted as developed by the previous 
authors, recognizing that some adjustments may improve the 
physical representation of the categorization along with its 
verification . 

B. PREDICTOR INTERPOLATION 

The interpolation method used to determine values of the 
predictors at the observation points is a natural bicubic 
spline curvilinear interpolation scheme. The locally de- 
veloped program, SPLIN, is available at the NPS W. R. Church 
Computer Center and was used throughout the study. 

C. REGRESSION SCHEME 

A stepwise multiple linear regression program, called 
BMDP2R (University of California, 1975) was chosen as the 
means for deriving the marine fog diagnostic equations. The 
polynomial regression scheme also availabe in the BMD com- 
puter program series was not used due to its inability to 
handle data sets of over 1000 cases. As each month has 
approximately 4000 synoptic observations and combined data 
sets had at least 7000 cases, the usefulness of the polynomial 
regression scheme applied to smaller segments of the data was 
unlikely to give accurate results. A non-linear, non-poly- 
nomial scheme also available in the BMD series uses excessive 
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computer memory for a 45-variable regression, and, therefore, 
was likewise not used. 

BMDP2R computes a sequence of multiple linear equations 

in a stepwise manner. At each step one variable is either 

added or removed from the previous step's equation as dependent 

on the F-to-enter and F-to-remove criteria. In the forward 

selection procedure the dependent predictor inserted in the 

equation is the one with the highest coefficient of partial 

2 

determination CAR ). At each step the regression procedure 
reevaluates the variables already in the equation, and_ may 
find that a variable important at an earlier stage may be 
less important at a later stage due to the intercorrelation 

of variables. In this procedure, before a variable is added, 
the variable already in the equation with the partial F 
value is dropped if this latter value is less than the maxi- 
mum F-to-remove value (.Wesolowsky , 1976). Selection proceeds 
until no candidate variables qualify or until all independent 
variables have been used. 

The BMDP2R program permits the specification of the 
tolerance level. If the tolerance is near zero, the variable 
being considered is close to being a linear combination of 
the variables already in the equation. Variables with low 
tolerances could cause computational difficulties in the 
method of calculation (Wesolowsky, 1976). Since tolerance is 
defined as one minus the coefficient of determination, the 
value of .01 was used as the minimum acceptable tolerance 
criterion (.University of California, 1975). 
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The approach taken in this study was to treat the y-inter- 
cept (Y-INTCP) as a variable and to use the BMDP2R program's 
default values of 4.0 and 3.9 as the values for the minimum 
F-to-enter and maximum F-to-remove criteria. 

D. VERIFICATION SCORING 

Two types of verification scores are used in the study to 
test the skill of different regression equations in describing 
the distribution of fog probabilities. The first type, 
exemplified by the Heidke Skill Score (HSS) and the Threat 
Score (TS), is used to test the effectiveness in specifying 
discrete occurrences of fog. The second type (Panof sky-Brier 
Probability Score) tests the accuracy of probabilities in 
estimating the likelihood of fog occurrence (PS) (Panofsky 
and Brier, 1958). 

The formulae used for the three scores (HSS , TS , and PS) 
are given in Appendix F. Currently the Threat Score is 
popularly used by the Techniques Development Laboratory, 
National Weather Service (Bermowitz and Best, 1978; Miller and 
Best, 1978), while the Heidke Skill Score has been traditional- 
ly used by meteorologists. The Panof sky-Brier Probability 
Score is suitable for evaluating the effectiveness of a 
predictor given in terms of a probability (0-100%) . 
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V . RESULTS 



A. REGRESSION EQUATIONS 

The sets of equations generated for the individual months 
of June, July and August 1976 and 1977 as well as combina- 
tions thereof, using categorization schemes FOGCAT I and 
FOGCAT II, are listed in Tables I -XII. They are shown in 
stepwise order and include Y-intercept , regression coeffi- 
cients and the amount of variance explained with the inclu- 

2 2 

sion or removal of each variable (R and AR ) . Only six 

steps are given in each table for purposes of brevity. In 

2 

all cases no two successive steps added a cumulative AR of 
greater than .005 after the sixth step until a variable was 
removed. Recognizing that variables which explain less than 
1% additional variance are normally not included in opera- 
tional regression equations, six steps were presented so 
that all variables having a possible effect on further devel- 
opment- could be identified. 

Since the climatological synoptic regime varies from 
month-to-month and from year-to-year the regression equa- 
tions vary accordingly. Several similarities are evident, 
however. First, in all cases for both FOGCAT I and FOGCAT II 
the first model output parameter to be entered is one of the 

heat flux terms, either EHF or SEHF . Second, either CLIMO 
or an interactive parameter containing CLIMO enters all the 
equations by step 3 with the exception of the equations 
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derived from August 1977 data and the equation for all 
months for FOGCAT I. Finally, the meridional component of 
the wind, v, or an interactive parameter containing v enters 
all of the equations no later than step 4 with the exception 
of the equations derived using FOGCAT I on the combined 
July 1976/77 data or FOGCAT I and II, on the August 1976 
and FOGCAT I on the August 1977 data. Additional similari- 
ties are evident at later steps but are not considered 
significant to this discussion. 

In all cases, the first three non-constant variables 
entered into the regression equations accounts for between 

71.9% and 92.5% of the summation of coefficients of partial 

2 

determination (AR ) determined by the entire set of model 
output predictors. 

As the magnitudes of the coefficients of determination 

2 

CR 1 vary greatly between those equations derived using the 
two different categorization schemes, a direct comparison 
of the effectiveness of the equations based on FOGCAT I 
versus those based on FOGCAT II is difficult. Although Table 
XIII indicates that the IAR^ for FOGCAT II is generally 
greater than this quantity for FOGCAT I, comparisons are more 
properly made using standard verification scoring methods, 
as in Section B following. 

The eight interactive parameters used in the regression 
analysis were formed by either multiplying the parameters 
found to be most significant by Quinn (1978) by climatology 
or by multiplying various of these most significant parameters 
by other significant parameters. The inclusion of these 
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interactive parameters allows climatology to enter the scheme 
at an earlier step. Originally, climatology was too highly 
correlated with the heat flux terms to enter the equations 
alone, and thus it was rejected by the tolerance criterion. 

Having identified the most probable meteorological indi- 
cators of fog occurrence by regression analysis, the feasi- 
bility of adding a persistence parameter was investigated. 

Due to lack of time, this investigation was limited to the 
July 1976 and 1977 data sets. The persistence parameter 
(PERS) was based only on those observations which had at 
least one synoptic report within a one degree latitude/longi- 
tude box of the reported position on the previous day. This 
requirement reduced the data set for July 1976 from 4391 to 
1368 reports and for July 1977 from 4198 to 1591 reports. 

Regression analysis, with PERS included, was then per- 
formed on the July data sets, both individually and in unison, 
using both categorization schemes. The sets of equations for 
these data sets are listed in Tables XIV-XVI . Persistence 
displayed a high correlation with EHE , SEHF and CLIMO and 
was rejected for inclusion in the equations by the tolerance 
criterion on this basis until the seventh step at the earliest. 
The similarities among equations which were noted earlier are 
decreased when PERS is included with only the heat flux term 
retaining its former importance in all cases. This instability 
introduced into the equations may be only the result of the 
reduced data set; the investigation of persistence as a fog 
indicator should not be abandoned until it is tested on a 
larger scale, both alone and in interactive combinations. 
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Fields of selected model output parameters, FTER and re- 



gression-computed fog probabilities for the arbitrarily 
selected observation time, 0000 GMT 3 August 1976, are shown 
in Figures 6-12. Considerable frontal activity in the mid- 
latitude central and western North Pacific Ocean is indicated 
by the sea-level pressure analysis (Figure 6). The two most 
important fog-related parameters, EHF and SEHF , are shown in 
Figures 7 and 8. Maximum positive values at the subtropical 
latitudes give way to a broad zonal band of negative and low 
positive values centered between 45-55N. Another important 
regression parameter, the algebraic value of meridional wind 
speed, v, is displayed in Figure 9. This field closely re- 
lates to the geostrophic/gradient wind implied by the sea- 
level pressure analysis. Figures 10 and 11 show the regres- 
sion-computed fog probabilities based on three variables for 

Fog Categorization Schemes I and II. The configuration of 
these fields resembles that of SEHF, the FNWC variable ex- 
plaining most of the variance in fog probability. Probabili- 
ties on the FOGCAT I version (Figure 10) appear to be about 
15% greater than those of FOGCAT II (Figure 11). The differ- 
ence is particularly noticeable in the subtropics where fog 
likelihood is climatologically small. The FTER probability 
field (Figure 12) shows a smaller range of values than the 
NPS probabilities, and is a closer fit to the FOGCAT II 
regression probabilities. 

B. VERIFICATION SCORING 

For purposes of verifying fog occurrence in an observation, 
major categories S and F were used in the FOGCAT I scheme, 
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and an assigned probability of fog occurrence of 100% was used 
in the FOGCAT IT scheme (.Appendices D and E)., An attempt was 
made to use 85% fog likelihood as also verifying fog occur- 
rence in FOGCAT II, but this criterion yielded significant- 
ly poorer results than 100% used alone, as had been found 
by Quinn (1978). This may indicate that some fine tuning of 
the assigned fog probabilities in FOGCAT II would yield 
better results. 

Optimum threshold probabilities were also computed for 
those equations based on the combined June data, combined 
July data, combined August data, and combined data for June, 
July and August 1976 and 1977. The threshold probability is 
that equation-specified probability which best separates the 
fog/no fog occurrences, hence yielding the best skill and 
threat scores (Bermowitz and Best, 1978). The threshold 
Heidke Skill Scores and Threat Scores presented in Tables 
XVII-XXI are those which maximize the indicated scores. The 
Panof sky-Brier Probability Score presented in these tables 
directly evaluates the accuracy of regression probabilities 
in estimating the likelihood of the event. 

The regression equations as derived by the BMDP2R program 
were not used in their entirety in the verification phase of 
this study nor were all the equations verified. Only the 
equations in the form containing the Y-intercept and the 
first three non-constant variables were used. It is possible 
to "overfit" the regression equation to the dependent data. 
Therefore only variables which bear a physical relationship 
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to the predictand should be included (Panofsky and Brier, 
1958). In all the equations verified, the first three vari- 
able predictors included had an apparently strong physical 
relationship with the occurrence of fog, while several equa- 
tions introduced terms with only a weak physical relation- 
ship to fog likelihood at step 5. 

The verification scores for the five data sets for which 
regression equations were evaluated and compared to the veri- 
fication scores for FTER and climatology are presented in 
Tables XV I I -XXI. Note that the scores for the regression 
equations as well as for FTER and climatology vary from the 
FOGCAT I to the FOGCAT II schemes. Since verification in 
each scheme is against the actual occurrence of fog (Category 
S or F in FOGCAT I or 100% in FOGCAT II), a relative compari- 
son of the effectiveness in categorizing fog occurrence 
between FOGCAT I and II for purposes of regression analysis 
is possible. This is not to imply that observations verify 
as fog/no fog equally with FOGCAT I and FOGCAT II; there are 
slight differences not considered important in this analysis. 

In four of the five cases verified, the skill and threat 
scores for the regression equations were higher for FOGCAT I 
than for FOGCAT II. In all verification scores, the FOGCAT 
II regression equation outperforms both climatology and FTER, 
except for the regression equation based on the two August 
data sets. The FOGCAT I regression equation outperforms 
climatology and FTER in skill and threat scores for each data 
set except the full June, July and August 1976 and 1977 set 
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where climatology is the top performer. In probability scor- 
ing, FTER outperforms the FOGCAT I regression equation in 
four of the five cases and climatology outperforms it in the 
other cases in the P-score. 

The variations of skill and threat scores for climatology 
from FOGCAT I to FOGCAT II may be taken as a measure of the 
differences in the verifying scheme (i.e. strong fogger, 
category S, and foggers, category F, for FOGCAT I and observa- 
tions assigned 100% probability in FOGCAT II). Thus, a mea- 
sure of the relative comparison of FOGCAT I and II regression 

equations might be obtained by noting the changes in the 
relative improvement over climatology. For example with 
reference to Table XVII, FOGCAT I skill/threat scores exceed 
climatologies skill/threat scores by .107/. 072, while for 
FOGCAT II the numbers are only . 082 /. 063 t indicat ing greater 
merit for FOGCAT I. Table XXII tabularizes this information. 
On this basis relative improvement over climatology is best 
using FOGCAT I for equations in Tables XVII, XVIII, XIX 
(except for threat score) and XXI, while FOGCAT II is better 
for equations in Table XX. 

When the persistence parameter is included in the regres- 
sion scheme for the two Julys, regression skill and threat 
scores generally show improvement (Tables XVIII-XIX). Con- 
sidering the reduced data set, at best, it is to be regarded 
as giving encouragement to the future refinement of a more 
comprehensive persistence parameter. 
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It is to be noted again that all regression equations 
apply to analysis time, hence are diagnostic. Further all 
are relative to dependent data. The data were considered 
too limited to expend even one month as an independent set 
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VI. CONCLUSIONS AND RECOMMENDATIONS 



From the evaluation of these diagnostic Model Output Sta- 
tistics regression equations, it can be seen that the MOS 
approach to prediction of marine fog holds great potential. 

The developed regression equations show improvement over both 
climatology and the FNWC advective fog forecasting product, 

FTER. However, it remains difficult to recommend, without 
qualifications, the use of either FOGCAT I or FOGCAT II 
schemes for the predictand. 

The comments that follow relate to the equation for each 
month (for two-year period); the equation for all six months 
combined are not considered to be operationally useful. 

Even though FOGCAT II equations score better than CLIMO and 
FTER (except for August) using all the scoring methods, 
threat, skill and probability, FOGCAT I equations show higher 
threat and skill scores than FOGCAT II and excel CLIMO and 
FTER in two of the scoring methods, threat and skill. It 
appears quite significant that FOGCAT I shows most skill 
relative to climatology and FTER as well and therefore this 
predictand categorization appears to be the best to use. 

Further development of FOGCAT I and II or another approach 
may well yield still further improvement in regression results. 

The following recommendations are offered as a guide to 
future work in this area: 



33 



1. In order to simplify the analysis procedure and the 
amount of data manipulation, parameters which have little or 
no effect on the outcome of the regression equation should 
be deleted. These include: SOLARAD , DDWW , PPW, PDW, SPW, 

SDW, WCP, THETAX, THETAR, STABX , ASTDX , ASTDR, CAPV , ADTSEA , 
AASTDX, AASTDR, and SSTA. See Appendix B. 

2. The development of a more realistic persistence para- 
meter should be continued and included in the regression 
scheme both alone and in interactive combinations. As a 
first approximation the fog probability assigned by the 
regression equation (without inclusion of a persistence para- 
meter) to each grid point of the previous day could be inter- 
polated to the actual observation positions of the present 
day. This would give a persistence value to each observation 
whether or not a report was available within a one degree 
latitude/longitude box on the previous day. These first 
approximations could then be weighted for those observations 
which had a report within the specified box on the previous 
day. In this manner the size of the data set would not be 
reduced as it was in this study, 

3. Additional data CP er haps one more year) should be 
added to further stabilize the regression analysis. 

4. A stratification of the set of regression equations 
by latitude, longitude, or by meteorological phenomena (i.e. 
wind direction, positive/negative thermal advection, or 
places where climatological fog frequencies are greater) is 
likely to improve accuracy. 
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5. The regression approach should be extended to regions 
other than the North Pacific Ocean and time periods other 
than summer. 

6. As FNWC predictive model output parameters fields 
become available for testing they should be incorporated into 
the regression equations for forecast intervals to at least 
48 hours . 
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TABLE XVI L Skill, 
data . 


threat and 


P-scores June 


1976 and 1977 


0 

A. FOGCAT I Categorization Scheme. R = . 


627. 




REGRESSION 

EQUATION 


FTER 


CLIMATOLOGY 


Skill Score 


* 

. 315 


.178 


. 208 


Threshold 


40 


14 


27 


Threat Score 


* 

. 338 


.258 


.266 


Threshold 


38 


21 


16 


P-Score 


.292 


.209 


* 

. 178 


Fog Probability (% 


) = 36.725 


- 2.095 (EHF) 






- .030 


(CLIMO) (SEHF) 


+ .226 (CLIMO) 



B. FOGCAT II 


Categorization Scheme 


. R 2 = 


.367. 




REGRESSION 

EQUATION 


FTER 


CLIMATOLOGY 


Skill Score 


.302* 


.207 


.220 


threshold 


28 


23 


27 


Threat Score 


* 

.322 


.250 


.259 


threshold 


26 


21 


22 


P-Score 


* 

.288 


.339 


. 300 



Fog Probability (%) = 31.154 - 2.956 (EHF) 

- .032 ( CLIMO ) ( SEHF ) + .023 (CLIM0)(V) 



* 

Indicates best score in each category 
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TABLE XVIII. Skill, threat and P-scores July 1976 and 1977 
data. 



o 

A. FOGCAT I Categorization Scheme. R = .670. 




REGRESSION 
EQUATION FTER 


CLIMATOLOGY 


Skill Score 


.384 .336 


. 303 


Threshold 


35 


35 


Threat Score 


.404* .355 


.360 


Threshold 


42 23 


28 


P-Score 


.330 .216* 


.242 


Fog Probability 


(%) = 34.898 - 1.663 (EHF) 






- .026 (CLIMO) (SEHF) + . 


.292 (CLIMO) 



B . FOGCAT 1 1 


Categorization Scheme 


. R 2 = 


.441 . 




REGRESSION 








EQUATION 


FTER 


CLIMATOLOGY 


Skill Score 


* 

.378 


.333 


. 302 


Threshold 


38 


35 


42 


Threat Score 


* 

.386 


. 346 


.347 


Threshold 


34 


23 


33 


P-Score 


* 

.299 


. 322 


. 322 



Fog Probability (%) = 37.131 - 4.037 ( EHF ) 

+ .029 OCL EMO ). ( V ) + .466 (EAIR) 

* 

Indicates best score in each category 
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TABLE XIX. Skill, threat and P-scores July 1976 and 1977 
data with persistence parameter. 



A. FOGCAT I Categorization Scheme. 


R 2 = .673. 




REGRESSION 






EQUATION 


FTER CLIMATOLOGY 


Skill Score .405 


. 334 


. 319 


Threshold 48 


28 


28 


* 

Threat Score .423 


. 371 


. 391 


Threshold 47 


25 


32 


P-Score . 356 


* 

.237 


.281 


Fog Probability (%) = 17.869 - 2.945 


(EHF) 




+ .510 (CLIMO) 


+ 1.113 (EX) 





B . FOGCAT I I 


Categorization Scheme 


. R 2 = 


.458. 




REGRESSION 








EQUATION 


FTER 


CLIMATOLOGY 


Skill Score 


* 

.394 


. 337 


. 320 


Threshold 


36 


33 


41 


Threat Score 


* 

.414 


. 365 


. 381 


Threshold 


36 


25 


33 


P-Score 


* 

. 323 


. 370 


. 350 



Fog Probability (%) = 28.167 - 3.956 (EHF) 

+ .642 (TX) - 130.026 (STABR) 



* 

Indicates best score in each category. 
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TABLE XX . Skill, threat and P-scores August 1976 and 
1977 data. 



A. FOGCAT I Categorization Scheme. 


R 2 = .660 




REGRESSION 

EQUATION 




FTER 


CLIMATOLOGY 


Skill Score .403* 




. 374 


. 146 


Threshold 45 




28 


26 


Threat Score .378 




. 345 


.219 


Threshold 40 




23 


14 


P-Score .279 




* 

. 175 


.202 


Fog Probability (%) = 38.781 - 


1.859 


( SEHF ) 




+ 1.059 


(V) - 


.708 (CAPU) 





B. FOGCAT II Categorization 


Scheme . 


R 2 = .420. 






REGRESS ION 








EQUATION 




FTER CLIMATOLOGY 


Skill Score 


* 

.402 




. 375 


. 136 


threshold 


46 




28 


25 


Threat Score 


. 368* 




. 341 


.208 


Threshold 


40 




24 


16 


P-Score 


.297 




.280* 


. 336 


Fog Probability 


(%) = 26.352 


- 2.105 


(SEHF) 






+ .727 


(V) - . 


067 ( SEHF ) ( V) 





* 

Indicates best score in each category. 
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TABLE XXI . Skill, threat and P-scores June, July and August 
1976 and 1977 data. 



2 

A. FOGCAT I Categorization Scheme. R = .646. 





REGRESSION 

EQUATION 


FTER 


CLIMATOLOGY 


Skill Score 


.214 


. 196 


* 

. 293 


Threshold 


27 


23 


44 


Threat Score 


.267 


.251 


* 

. 319 


Threshold 


16 


21 


41 


P-Score 


.304 


* 

.201 


.208 


Fog Probability 


(%) = 31.966 - 1 


.559 ( SEHF) 






+ .561 (V) 


+ .215 (CLIMO) 





B. FOGCAT II Categorization 


Scheme. R^ = .400. 




REGRESSION 

EQUATION 


FTER CLIMATOLOGY 


Skill Score 


* 

. 370 


.321 .294 


Threshold 


32 


34 40 


Threat Score 


* 

.382 


.335 .337 


Threshold 


31 


23 32 


P-Score 


* 

.283 


.322 .318 


Fog Probability 


(%) = 30.406 
+ .029 


- 2.811 (SEHF) 

(CLIMO) (V) + 1.530 (SHF) 


* 

Indicates best 


score in each 


category . 
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TABLE XXII: Comparison of Skill and Threat Scores for 

the Regression Equations and Climatology 



Period 


FOGCAT I 


FOGCAT II 


Skill Score 
Difference 

(Regression 
Equation - 
Climatology) 


Threat Score 
Dif f er ence 
( Regress ion 
Equation - 
Climatology) 


Skill Score 
Difference 

(Regression 
Equation - 
Climatology) 


Threat Score 
Difference 
( Regression 
Equation - 
Climatology) 


June 
76 & 77 


.107 


.072 


.082 


.063 


July 
76 & 77 


.081 


.044 


.076 


.039 


July 76 
& 77 with 
persis- 
tence 


.086 


.032 


.074 


. 033 


August 
76 & 77 


. 257 


. 159 


.266 


. 160 
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* 
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Figure 1. Study area on polar stereographic projection 




Figure 2. Fleet Numerical Weather Central's 63x63 

grid, with outline of North Pacific Ocean 
rectangular grid area used in study . See 
text . 
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Figure 3. Percent frequency of occurrence of fog in any form 
June (Guttman, 1978). 
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Figure 4. Percent frequency of occurrence of fog in any form 
July (Guttman, 1978). 
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Figure 5. Percent frequency of occurrence of fog in any form 
August (Guttman, 1978). 
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, 



Figure 6. National Meteorological Center's Sea-Level Pressure Analysis, 0000 GMT 
3 August 1976, North Pacific Ocean. 




, 
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Figure 8. FNWC's Sensible/Evaporative Heat Flux, SEHF , 0000 GMT 3 August 
North Pacific Ocean. 
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Figure 10. FNWC ' s Advective Fog Probability FTER, 0000 GMT 3 August 1976 
North Pacific Ocean. 



VS-/0, -? I 




, * fr 



68 




X 



X 






CD 

t> 1 

CD > 

rH 

t> 

-M CM 
0} t> 

3 * - 

bflO C 

$ Cj 

< + 0 
O 

a p o 

o x 

U P o 

Hh C/3 -H 
«H 

^ lO -H 

0 O O 

> rH CS 

■H • P 

(H CM 

0 i x: 

Q CM -M 

m ^ 

C CO 

O 

*H CD 

-P CM 

3 II 

C? 

P >> rH 

-P 

£ -H 4J 
OH Cfl 
■ri *H 3 
Cfi P b£) 
a: d p 
0 P <3 
Sh O 
&£) Sh CO 
0 p 
p — H 



* 



O 



CD 

05 






CO HH 

P HH 

S <5 

O CJ 
in a 
O 
P 

+J - 
■H CCS 
rH +J 

•h as 

p a 

ccS 

X2 t> 

° £ 

Sh 05 

p rH 

aflTO 
O fl 
p aS 



M 



a 



o 

o 

o 

o 



p 

3 

CO 

X 

> 



CD 

O 



■r 



0 

u 

3 

cJ D 

•H 

P 



4 

o<- 






69 



APPENDIX A 



LAND STATIONS USED IN STUDY 
(U.S. Air Force, 1972) 



Station 

Number 


Latitude 


Longitude 


Surface 

Elevation 


32174 


45.2 


N 


147.9 


E 


38 


32186 


46.2 


N 


150.5 


E 


73 


32195 


46.9 


N 


151.9 


E 


26 


32207 


43.3 


N 


153.3 


E 


55 


32213 


50.9 


N 


156.7 


E 


42 


32215* 


50.7 


N 


156.2 


E 


unknown 


32217 


50.0 


N 


155 . 4 


E 


11 


32559 


53.1 


N 


160.0 


E 


88 


32594 


51.5 


N 


156.5 


E 


6 


32618 


54.9 


** 

N 


166.2 


** 

E 


19 


70454 


51.9 


N 


176.7 


W 


14 



Station not listed in reference. Lat/long values 
listed are as given with each synoptic report 
furnished by FNWC. 

Location given slightly removed from actual loca- 
tion of 55.2 N latitude, 166.0 E longitude to 
accommodate study area. 
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APPENDIX B 



OUTPUT PARAMETER DESCRIPTIONS 



Source Model 

Name 

Symbo 1 Description Units 

A . Analysis Parameters (FNWC's Mass Structure Model) 



PS Sea Level Pressure (mb) 

Analysis of observed sea level 
parameter. 

TAIR Surface Air Temperature (°C) 

Analysis of observation-level 
air-temperature. 

EAIR Surface Vapor Pressure (mb) 

Analysis of observation-level 
vapor pressure derived from the 
dew point. 

T925 925 mb Air Temperature (°C) 

Analysis of 925 mb air tempera- 
ture . 

TSEA Sea Surface Temperature (°C) 

Once-daily analysis of observed 
sea-surface temperature. 



B. P.E. Parameters (FNWC's Primitive Equation Model) 

TX Surface Air Temperature (°C) 

Derived from surface air and 
potential temperatures, boundary 
layer depth, upper-level winds 
extrapolated to surface, air 
density, drag coefficient, gusti- 
ness factor, and empirical con- 
stants . 



EX Surface Vapor Pressure 

Derived from model's mixing ratio. 



(mb) 



SOLARAD Solar Radiation 

Calculated absorption of incoming 
short-wave (solar) radiation 
(positive downward) 



(gm-cal/ 

2 u \ 
cm hr ) 
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EHF 


Evaporative Heat Flux 
Derived using air density, drag 
coefficient, extrapolated winds, 
and mixing ratios. 


(gm-cal/ 
cm^ hr) 


SERF 


Sensible Plus Evaporative Heat Flux 
SEHF = SHF + EHF 


( gm-cal/ 
cm hr) 


SHF 


Sensible Heat Flux 
Recovered from SHF = SEHF - EHF. 
Originally derived by FNWC using 
drag coefficient, extrapolated 
winds, surface air temperature, TX, 
density, and constants. 


( gm-cal/ 
cm hr) 


THF 


Total Heat Flux 

THF = SEIIF - SOLARAD + LW, where 
LW is the heating due to long-wave 
(terrestrial) radiation. 


(gm-cal/ 
cm hr ) 


Marine 


Wind Model (FNWC) 




WWW 


Marine Wind Speed 


(knots ) 


DDWW 


Marine Wind Direction 
Both variables derived from a 
dynamic balancing of surface wind 
and sea-level pressure. 


( degrees/ 
10) 


Spectral Ocean Wave Model (S .O.W.M . ) (FNWC) 




HW 


Significant Wave Height 


( feet ) 


PPW 


Primary Wave Period 


( sec ) 


PDW 

SPW 


Primary Wave Direction 
Secondary Wave Period 


( degrees/ 
10) 

( sec) 


SDW 

WCP 


Secondary Wave Direction 
Probability of White Caps 


( degrees/ 
10) 

(percent ) 


Other 


Model Output Parameters (FNWC) 




SSTA 


Sea Surface Temperature Anomaly 
Calculated anomaly of sea-surface 
temperature from the mean of the 
day as interpolated from the 
monthly mean values. 


(°C) 



F. Derived Parameters 



u 


Zonal Wind Component 
u = -WWW sin(DDWW • 10) 


(m/sec) 


V 


Meridional Wind Component 
v = -WWW cos(DDWW • 10) 


(m/sec ) 


CAPU 


I Directional Wind Component 

CAPU = -u • sin(LNGA) - v • cos(LNGA), 

(Haltiner, 1971). 


(m/ sec) 


CAPV 


J Directional Wind Component 

CAPV = u ♦ cos(LNGA) - v • sin(LNGA) , 

(Haltiner, 1971), where LNGA = 

-10 - (I,J point longitude). 


(m/sec ) 


THETAX 


Potential Temperature X 
Derived using PS, TX. 


(°K) 


THETAR 


Potential Temperature R 
Derived using PS, TAIR. 


( °K ) 


STABX 


Stability X 

Derived using (THETAX- ( THETA of T925); ] 
/(PS-925). Value greater than zero 
indicates absolute instability. 


( °K/mb ) 


STABR 


Stability R 

Derived using (THETAR- (THETA of T925) ] 
/(PS-925). Same value effect as STABX. 


( °K/ mb ) 


ASTDX 


Air-Sea Temperature Difference X 
ASTDX = TX - TSEA 


(°C) 


ASTDR 


Air-Sea Temperature Difference R 
ASTDR = TAIR - TSEA. 


(°C) 


ADTSEA 


Advection of TSEA 
Formulae and notes below. 


(°c/ 

Hour) 


ADTX 


Advection of TX 
Formulae and notes below. 


(°c/ 

Hour) 


ADTAIR 


Advection of TAIR. 
Formulae and notes below. 


(°c/ 

Hour) 


AASTDX 


Advection of ASTDX 
Formulae and notes below. 


(°c/ 

Hour) 
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AASTDR 


Advection of ASTDR 
Formulae and notes 


below. 


C°c/ 

Hour) 


Climatological Parameters 






CLIMO 


National Climatic Center Fog 
Frequency Climatology 


(percent/ 

100) 


PERS 


Persistence of Fog 


Likelihood 


( percent / 
100) 


Interactive Parameters 






CLISEHF 


CLIMO * SEHF 




(% gm-cal/ 
q 

100 cm hr) 


CLIV 

VSEHF 


CLIMO * v 
v • SEHF 




(m %/ 

100 sec) 

0-00 gm-cal/ 
cm sec hr) 


SEHFSQ 


SEHF • SEHF 




(gm -cal / 

n p 

cm" hr) 


CLIEHF 


CLIMO * EHF 




(% gm-cal/ 
100 cm^hr) 


CLIADT 


CLIMO • ADTAIR 




(% °c/ 

100 hr) 


EHFSQ 


EHF • EHF 




(gm -cal / 

2 2 
cm hr) 


EHFADT 


EHF • ADTAIR 




(gm-cal °C/ 

2 , 2v 

cm hr ) 



Advection Formulae and Conditions: 

For the advection of a quantity (R) the formula, ADQ = 
-V*V(Q), was used in the following finite difference form: 



ADC ! - - Tjr [CAPD -<‘5l + l- Q I-l>J + CAPV -<Qj + l-‘Wl 1 ' 

where RMAP = (l+sin(60) )/ (l+sin( latitude) ) 

and DM = [ (2) • (6 . 37- 10 6 ) • ( l+sin(60) ) ] /31 . 205 

(31.205 = grid mesh lengths, pole to equator, on FNWC’s 
I,J grid). 

In the temperature advection calculation for point C, using 
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the five grid points illustrated below, one or two of the 
points, namely, A, B, D, or E , may be outside the study area. 
In the bogusing method suggested by Mr. Leo Clarke, FNWC, 
when a non-center point (e.g., point A) was judged to prob- 
ably produce a land and/or dissimilar sea area influence on 
the resulting advection, the center point value (point C) 
was substituted for it. This "bogusing" is necessary to 
maintain a "purely" marine characteristic in the resultant 
parameter value. 

The set of study area boundary grid points used for 
bogusing (and some double bogusing) are depicted in Fig. 1. 
The upper case letters denote the points whose values were 
used for bogusing. The lower case letters mark the posi- 
tions of the adjacent points being bogused. The only bound- 
ary point close enough to land to give concern is the one 
near station number 32594 (south tip, Kamchatka Peninsula). 
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APPENDIX C 



Abridged version of internationally used weather code 
figures and definitions for reporting present and past 
weather and low clouds in the surface synoptic report 
(U. S. Departments of Commerce, Defense, and Transporta- 
tion , 1969 ) . 

Present Weather 



Code 

Value 


Definition 


00-03 


Characteristic change of the state of the 
sky (cloud) during the past hour. 


04-09 


Haze, dust, sand, or smoke. 


10 


Deep light fog. 


11-12 


Shallow heavy fog. 


13-17 


Lightning, thunder, or precipitation within 
sight, not reaching the ground. 


18-19 


Squall(s), funnel cloud(s) during the past 
hour . 


20 


Drizzle during the past hour. 


21-23 


Rain, snow, or rain and snow during the past 
hour . 


24 


Freezing drizzle during the past hour. 


25-27 


Shower(s) during the preceding hour. 


28 


Fog during the past hour. 


29 


Thunderstorm during the past hour. 


30-39 


Duststorm, sandstorm, drifting or blowing snow 


40 


Fog at distance, but not at station, during 
the past hour (visibility less than 1 km). 


41-49 


Deep heavy fog at the time of observation 
(visibility less than 1 km). 


50-59 


Drizzle, or drizzle and rain. 


60-63 


Slight to moderate rain. 


64-65 


Heavy rain. 


66 


Slight freezing rain. 


67 


Moderate or heavy freezing rain. 


68 


Slight rain or drizzle and snow. 


69 


Moderate or heavy rain or drizzle and snow. 
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70-79 Solid precipitation not in showers. 

80-99 Showery precipitation or precipitation with 
current or recent thunderstorms. 



Past Weather 



0 Cloud covering i or less of sky throughout the 
period. 

1 Cloud covering more than i of sky during part of 
the period. 

2 Cloud covering more than i of sky throughout period. 

3 Sandstorm, or duststorm, or blowing snow. 

4 Heavy fog, thick haze, or smoke. 

5 Drizzle 

6 Rain 

7 Snow, rain and snow. 

8 Shower (s). 

9 Thunderstorm, with or without precipitation. 

Low Cloud Type 

0 No low clouds. 

1 Ragged cumulus of fair weather. 

2 Generally towering cumulus. 

3 Cumulonimbus without cirriform or anvil tops. 

4 Stratocumulus formed by cumulus spreading out . 

5 Stratocumulus not formed by cumulus spreading. 

6 ■ Stratus or f ractostratus . 

7 Fractostratus of bad weather. 

3 Cumulus and stratocumulus, with bases at different 
levels . 

9 Cumulonimbus with cirriform top. 

/ Low cloud obscured. 

Visibility 



90 Less than 50 m 

91 0-199 m 

92 200-499 m 

93 500 m - 0.99 km 

94 1 - 1.99 km 

95 2 - 3.99 km 

96 4 - 9.99 km 

97-99 equal to or greater than 10 km 
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APPENDIX D 



FOGCAT I Categorization Scheme 
(Van Orman and Renard, 1977) 

1. Groupings and symbols used in FOGCAT I categorization 
scheme . 



Present 


Weather (ww) 


Past Weather (W) 


Low 


Cloud (CL) 


Symbol 


Associated 


Symbol 


Associated 


Symbol Associated 




ww Codes 




W Codes 




CL Codes 


41G 


41-49 


4 


4 


6 


6 


10G 


10,28,40 


4,5 


4,5 


5,7 


5,7 


11G 


11,12,20,24 


5G 


0, 1,2,5 


B 


/ 


50G 


50-59 


2G 


0,1,2 


* 


any CL not 
listed above 


60G 


60-63,66,68 


* 


any W not 
listed above 






* 


any ww not 
listed above 











G = Group B = Low clouds obscured 

2. Scheme for categorizing observations according to like- 
lihood of fog (FOGCAT I). 



Major 


Sub- 


Present 


Past 


Low Cloud 


Category 


Category 


Weather 

(ww) 


Weather 

(W) 


Type 

(CL) 


Strong 


al 


41G 


4,5 


6 


Foggers 


a2 


T T 


H 


B 


= S 


bl 


Tt 


2G 


6 




b2 


IT 


TT 


B 




b3 


10G 


4,5 


6 




b4 


IT 


TT 


B 




cl 


41G 


* 


6 




c2 


IT 


4,5 


5,7 




c3 


1 T 


TT 


* 




c4 


10G 


2G 


6 



Foggers 


dl 


41G 


* 


B 


= F 


d2 


TT 


2G 


5,7 




d3 


50G 


4 


6 




d4 


TT 


u 


B 




d5 


10G 


* 


6 




d6 


TT 


2G 


B 




d7 


TT 


4 


5,7 
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Past/Weak 
Foggers 
= P 



Maybe 
Foggers 
= M 



el 


41G 


* 


5,7 


e2 


1! 


2G 


* 


e3 


50G 


4 


5,7 


e4 


10G 


* 


B 


e5 


TT 


4 


* 


e6 


11G 


4 


6 


el 


II 


II 


B 


f 1 


41G 


* 


* 


f 2 


50G 


4 


* 


f 3 


10G 


5G 


5,7 


f 4 


11G 


1 1 


6 


f 5 


II 


II 


B 


gl 


10G 


* 


5,7 


g2 


ii 


5G 


* 


g3 


11G 


* 


6 


g4 


II 


4 


5,7 


g5 


f! 


II 


* 


hi 


M 


* 


B 


h2 


If 


5G 


5,7 


il 


10G 


* 


* 


i2 


11G 


* 


5,7 


jl 


IT 


5G 


* 


kl 


50G 


5G 


6 


k2 


GOG 


4 


6 


11 


IT 


II 


B 


12 


* 


1 1 


6 


ml 


50G 


* 


6 


m2 


* 


4 


B 


nl 


* 


5G 


6 


ol 


60G 


4 


5,7 


o2 


* 


4 


5,7 


Pi 


60G 


4 


* 


p2 


* 


* 


6 


p3 


* 


4 


* 


ql 


50G 


5G 


B 


q2 


n 


II 


5,7 


q3 


ii 


* 


5,7 


q4 


it 


5G 


* 


q5 


11G 


* 


* 


q6 


60G 


5G 


6 


q7 


ii 


* 


6 


rl 


50G 


* 


B 


r2 


n 


* 


* 


r3 


60G 


5G 


B 


r4 


it 


ii 


5,7 


r5 


ii 


* 


B 


r6 


ii 


* 


5,7 
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Non- 


ul 


60G 


5G 


* 


Foggers 


u2 


n 


* 


* 


= N 


vl 


* 


5G 


B 




v/1 


* 


* 


B 




xl 


* 


* 


5,7 




yi 


* 


5G 


5,7 




y2 


* 


5G 


* 




zl 


* 


* 


* 


3. Major 


fog 


categories 


and sub-category 


fog groups with 




associated Fog 


Probabilities as 


a function of 



visibility . 







Visibility 


Code Values 






96-99 


94-95 


90-93 


Fog 




(Good) 


(Fair ) 


(Poor) 


Prob abi! 


Major 

Category 




Subcategory 


Groups 










a 


100.0 






a 


b 


96.6 




a 


b 


c 


93.1 


S 


b 


c 


d 


89.7 




c 


d 


e 


86.2 




d 


e 


f 


82.8 




e 


f 


g 


79.3 




f 


g 


h 


75.9 


F 


g 


h 


i 


72.4 




h 


i 


j 


69.0 




i 




k 


65.5 

62.1 




“k” 


k 


£ 


58.6 




£ 


£ 


ra 


55.2 




m 


ra 


n 


51 . 7 


P 


n 


n 


o 


48.3 




o 


o 




44.8 




P 


P 




41.4 




q 


q 


q 


37.9 


M 


r 


r 


r 


34.5 


not 


s 


s 


s 


31.0 


used 


t 


t 


t 


27.6 








u 


24.1 




u 


u 


V 


20.7 




V 


V 


w 


17.2 




w 


w 


X 


13.8 


N 


X 


X 


y 


10. 3 






y 


z 


6.9 




y 


z 




3.4 




z 






0.0 
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APPENDIX E 



FOGCAT II Categorization Scheme (Quinn, 1978) 



Fog 

Present Weather Past Weather Visibility Probability 

(w) cwj (vis) m 



0-9 


- 


- 


0 


10-12 


— 


- 


100 


13-19 


— 


— 


0 


20 


o 

1 

CD 

* 

CO 

1 

o 


— 


0 


20 


4,5 


90-95 


35 


20 


4,5 


96-99 


0 


21-24 


0-3, 6-9 


- 


15 


21-24 


4,5 


90-95 


60 


21-24 


4,5 


96-99 


15 


25-27 


0-3, 6-9 


- 


0 


25-27 


4,5 


90-95 


10 


25-27 


4,5 


96-99 


0 


28 




- 


100 


29-39 


— 


— 


0 


40-49 


— 


— 


100 


50-59 


0-3, 6-9 


- 


35 


50-59 


4,5 


90-95 


85 


50-59 


4,5 


96-99 


35 


60-69 


0-3, 6-9 


- 


15 


60-69 


4,5 


90-95 


60 


60-69 


4,5 


96-99 


15 


70-79 


— 


— 


0 


80-89 


O 

1 

CD 

CO 

1 

o 


— 


0 


80-39 


4,5 


90-95 


10 


80-39 


4,5 


96-99 


0 


90-99 




— 


0 



Dash indicates that the particular category was not con- 
sidered in the assignment of a fog probability. 
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APPENDIX F 



Verification Score Formulae (Quinn, 1978) 



EVENT ESTIMATED 



YES 

YES A 

EVENT 

OBSERVED 

NO B 

Heidke _ FC - EX 

Skill Score T - EX 



NO 

C Total (T) = A+B+C+D 

No. of Correct 
D Forecast FC = A+D 

Range: I 2 — 9 < HSS < 1 

a + C 



where EX = ( A+BHA+C ) + _ ( D+B ) ( D+C ) ^ No> of expectec i 

correct forecasts due to chance. 

Threat Score = Range: 0 £ TS _< 1 

Both scores indicate more skill with larger 
positive values. 

The Probability Score (PS) is from that given by 
Panofsky and Brier (1958) and may be written as (Renard, 
1975): 



P-Score 



where 





n 

o 


n l 




2 

N 


i£ 


pi 2+ E 


a-Pj> 2 ] 




i=l 


j=i 








Range : 



N = Total number of cases 
n Q = Total number of non-events 

Pp = Associated probability value for the non-event 
n^ = Total number of events 

Pj = Associated probability value for the event 
The closer to zero, the greater the skill. 
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